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We present an example of a limit cycle, i.e., a recurrent flow-line of the beta- function vector 
field, in a unitary four-dimensional gauge theory. We thus prove that beta functions of four- 
dimensional gauge theories do not produce gradient flows. The limit cycle is established in 
perturbation theory with a three-loop calculation which we describe in detail. 
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1. Introduction 



A necessary prerequisite for the complete understanding of quantum field theory (QFT) 
is the appreciation of its possible phases. In some cases a phase may be out of direct 
computational reach, e.g., the confining phase of QCD, while in others one may be able 
to use perturbation theory to gain an understanding of the dynamics of the theory. For a 
long time the only perturbatively accessible phase of QFT has been presumed to be that 
of a theory at a conformal fixed point, where, e.g., correlators exhibit power-law scaling. 

Recently, the existence of renormalization-group (RG) limit cycles was established by 
us in d = 4 — e spacetime dimensions with a three-loop calculation in a unitary theory of 
scalars and fermions (l}]3). Theories in d = 4 — e are of course unphysical, but working 
with them has always been useful in the study of properties of the RG [4], in the sense 
that RG effects found in such theories have invariably been shown to have counterparts in 
more physical cases. It was therefore suggested by our results that limit cycles should also 
occur in integer spacetime dimensions. In the present note we show that this is indeed the 
case in a four-dimensional unitary gauge theory. 

This new feature of the RG gives rise to an obvious question: "what phase of QFT is 
described by a limit cycle?" It follows from the work of Jack and Osborn |5| that theories 
that live on limit cycles may be CFTs. As we show in [6] this is indeed the case for the 
limit cycle we present below. Thus, although beta functions admit limit cycles, theories 
that live on these cycles are conformal. 

The existence of recurrent trajectories in the RG has implications for the c-theorem. 
This theorem reflects the intuition that coarse-graining reduces the number of massless 
degrees of freedom of a QFT, and it comes in different versions, as explained, e.g., in [7]. 
The strong version, i.e., that there exists a scalar function of the couplings c, along any RG 
flow, that obeys dc/dt < 0, with t the RG time and the inequality saturated only at fixed 
points, was proved long ago for QFTs in d = 2 [8] , and has been elaborated on heavily in 
the literature. Soon thereafter it was suggested that a strong c-theorem should be true in 
d = 4 as well [9], and that was indeed shown to be the case at weak coupling [5,10 , at least 
when renormalization effects of certain composite operators are not of relevance [3 11 . A 
proof of the four- dimensional version of the weak version of the c-theorem was recently 
claimed [12] (see also [13]), i.e., that there is a c-function such that if two four-dimensional 
CFTs are connected by an RG flow, then cuv > c ir- Similar ideas were used in an attempt 
for a proof of the weak version of the c-theorem in d = 6 
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We hasten to remark that the existence of limit cycles in the beta-function vector field 



1 



does not contradict intuition derived from the c-theorem. In particular, the quantity c that 
satisfies a c-theorem is constant even on limit cycles, and is expected to have the same 
monotonic behavior when it flows from a UV fixed point or limit cycle to an IR fixed 
point or limit cycle. However, the existence of RG limit cycles obviously demonstrates that 
beta-function flows are not gradient flows. 

The outline of the paper is as follows. In the next section we present our example. 
We describe in detail the three-loop calculation that establishes the limit cycle, and we 
show that the dilatation current of the theory on the limit cycle is well-defined and free of 
anomalies. In the last section we conclude and mention a few open questions. 



2. The 4d example 

In this section we describe in detail the first example of a limit cycle in d = 4. 
2.1. The theory 

Our theory has an SU(3) gauge group with two singlet real scalars, 4>i an d 02, two pairs 
of fundamental and antifundamental active Weyl fermions, (^1,2, ^1,2), as well as §(29 — 3e) 
pairs of fundamental and antifundamental sterile Weyl fermions. The kinetic terms are 
canonical and the interactions are given by[j] 

*) 3 (JH* *) (; 

In contrast with the active Weyl spinors, the sterile ones do not interact with the scalars, 
but they do interact with the gluons through their kinetic terms. One needs sterile fermions 
in order to get a perturbative fixed point for the gauge coupling, a la Banks-Zaks [2 15 . 
The smallest value of e for which our theory is physical is e = |, but we will treat e as an 
expansion parameter and take e — > | at the end. As we will see, our perturbative results 
can be trusted in this limit. 



1 The beta functions for all couplings in this theory can be found at http://het.ucsd.edu/misc/4D_ 



betas2s!2f .m 
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The most general virial current in our theory i^] 



(2.1) 



P ir For 



where Q a b is antisymmetric and Pij anti-Hermitian, i.e., Qb a = —Qab and = 
compactness we have denoted by ip 34: the two antifundamentals ipi 2 . By gauge invariance 



Pji = for % = 1,2 and j = 3,4. All these constraints are satisfied by 



Q 




and 



P 
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(2.2) 



-p 7 + ip 8 ipi J 



The virial current (2.1) contains a fermionic part, something that can lead to an ABJ- 



like anomaly 16 for the dilatation current. In this case, the virial current would have an 



extra contribution to its anomalous dimension^] beyond the one calculated from its two- 
point function. This is not allowed by conformal invariance |6j. Therefore, we expect that 
a limit cycle solution should have the property that the virial current be not anomalous. 
Consequently, a consistent limit cycle in a gauge theory with an SU(n > 3) gauge group 
and fundamental and antifundamental fermions should have 

Tr P = 0. 

This condition provides a powerful check on our computations. 
2.2. The three-loop computation 

It is convenient to rewrite compactly the interactions in V as 

v = ^Kbcd4>a4>b4>c4>d + {\ya\ijM%^j + h.c). 

Here, again, we are using the compact notation for the Weyl spinors, with 4 standing 
for ^i )2 . To find a limit cycle we must exhibit solutions to 

P a (g,y,\) = 0, 

Pa\ij(9, A) = —Qa'aVa'lij ~ Pi'iVa\i'j ~ Pj'jValij' > (2-3) 
Pabcdid, V, A) = —Qa'a^a'bcd ~ Qb'b^ab'cd ~ Qc'c^abc'd ~ Qd'd^abcd' > 



2 Lower case indices from the beginning of the roman alphabet are indices in flavor space for scalar fields, 
while lower case indices from the middle are indices in flavor and gauge space for Weyl spinors. 
3 As is also the case, e.g., for the axial current 17 Appendix C]. 
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that do not require zero {3 a \ij and/or fi a bcd- This requires both determining the values of 
the coupling constants and of the matrices Q and P for which the equations are satisfied. 
It would appear, naively, that the system of equations (2.3) has more unknowns than equa- 
tions, due to the presence of the unknowns Q a b and Py, and is thus ill-defined. However, in 
searching for particular solutions, one is free to set some coupling constants to zero. This 
is accomplished by using the freedom to redefine the scalar fields by an 0(2) transforma- 
tion and the active Weyl spinors by a U(2) x U (2) transformation, with the concomitant 
redefinition of coupling constants. Note that a coupling may become zero without its beta 
function becoming zero, since the couplings are not exclusively multiplicatively renormal- 



ized. Hence, the number of unknowns in (2.3) is reduced and we obtain a well-defined 
system with equal numbers of equations and unknowns. 

As in Ref. [3] we can calculate the entries of Q and P on a limit cycle in an expansion 
in e. To that end, we expand in the small parameter e the couplings, 

9 - / jT £ 2 > Va\ij - y a \ij £ 2 > A <>bcd ~ A abcd £ > 

n>l n>l n>l 

and the unknown parameters in the virial current, 

Qab = ^ Q { ab £n i P ij = ^2 P if £n ' 

n>3 n>3 



and we solve Eqs. (2.3) order by order in e. The lowest order entries in Q and P are of 
order e 3 , for at lower orders in e, corresponding to one- and two-loop orders in perturbation 
theory, the beta functions produce a gradient flow [H]. 

To establish a limit cycle we have to compute the e 3 -order terms in the e expansion of 
the parameters of the virial current. For a complete calculation of these we need the two- 
loop beta function for the quartic coupling, the three-loop beta function for the Yukawa 
coupling, and the four-loop beta function for the gauge coupling. To see why, let us explain 
how the e expansion works. 

The e expansion of the beta functions can be written schematically as 



g 2 

a n>l 



+ f^\g {1 \y {1 \x {1 \g {2 \y {2 \x i2) -,g {3 \y {3) ) £7/2 + --- 

+ /f (g {l) J l \ \ {l \g {2 \y {2 \ X {2) ; g i3 \y {3) ) £7/2 + ■■■ 
P x = E lf }£n+1 = fx'i^J 1 '; A (1) )e 2 + f?(g w ,y m , \ m ,g {2 \y®; X (2) )e 3 + ■■■ . 



n>l 
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Note that the gauge-coupling beta function is divided by g 2 . This way systems of equa- 
tions obtained at a specific e order contain the same coefficients in the e expansion of the 
couplings and can thus be solved simultaneously. All couplings, /'s, and beta functions 
carry flavor indices which we omit for brevity. It is important to realize that both the one- 
and the two-loop order of (5 9 contribute to fg 1 ^, for we are fixing the gauge coupling to 
a point a la Banks-Zaks. The first step is to simultaneously solve fg = and fy = 0, 
a system of nonlinear equations from which we get a set of solutions {(g^\ y^)} ■ Each 
solution in this set is then used to solve fj^' = 0, another system of nonlinear equations, 
which also gives a set of solutions {A 1 - 1 )}. At this point we can discard solutions with 
complex A^'s — those correspond to nonunitary theories — and construct the set of solu- 
tions iS = {(g^\ A^ 1 - 1 )}. The determination of the unknowns in /i n ~ 2 ' ) requires solving 
simultaneous linear equations, and so we have a unique solution for each element of S. 

(2) (2) 

At the next step we use solutions in S to simultaneously solve fg — and f y = for 
the unknowns g^ and y^ 2 \ which are thus uniquely determined. These are used in 
from which is determined, and then we consider fg 3 ^ and fy . These two functions 
receive contributions from the (n < 4)-loop orders of ft 9 and the (n < 3)-loop orders of 
(3 y . At this level we also have to take Q and P into account, i.e., we have to see if there 

(3) 

are solutions in the set S that can lead to solutions of the linear equations fg — and 
f!y°> = Q y + p y with Q and/or P nonzero. An indication of which solutions in S may lead 
to non- vanishing Q or P is that, already at the previous order, the beta functions for the 
coupling constants that were set to zero do not vanish. 

Now, the two-loop Yukawa and scalar coupling beta functions and the three-loop gauge 



beta function can be found in the literature 18 . To establish the non-vanishing of Q 
or P at lowest order, e 3 , the three-loop Yukawa beta function and the four- loop gauge 
beta function are required. Fortunately, a complete calculation of these beta functions is 
not needed. We parametrize the beta functions at these orders by summing all possible 
monomials of coupling constants of appropriate order with arbitrary coefficients c n . Then, 

(3) (3) 

by solving the set of linear equations fg — and fy = Qy + Py, we determine which 
of these coefficients are involved in the determination of Q and P. There is a one-to- 
one correspondence between each monomial in the beta functions and a three- or four- 
loop Feynman diagram. Thus, rather than computing some 1200 three-loop diagrams for 
the Yukawa beta function and a larger number of four-loop diagrams for the gauge beta 
function, we find that only a small number of diagrams needs to be computed. 

For the present model, following the procedure outlined in the previous paragraph, we 
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find, to lowest order in e, that the point 

_ 219a/ 92 534(11 430 301-30 212v / T9^5) 3/2 , ■ 24 v / 3559 fx 
y i ~ 27828 258 757 7T£ ~*~ 1 3559 71 v e + ' ' ' 



_ 8^17795 /7 | _ 16^10677 /7 , 

2/ 4 - 3559 e + , V5 ~ 3559 ^V^ + 



\ _ -3(4 177004+11 781 V19 370) 2 , \ _ -75(93 964+1245\/19 370) 2 , /- 

Ai- 7819123 7T £+•--, A 2 - 7819123 7T £ + ■ • ■ , (2.4) 

\ _ 1743(9y / 19370-676) 2 , \ _ -249a/ 78(11 430 301-30 212^19 370) 2 , 

3 ~~ 7819123 7T £ -h • • • , A 4 — 7819123 7T £ "h • • • , 



x _ -63^78(11 430 301-30 212v / l9370) _ 2 ,_ , „ _ 6x778298 ^ fZ , 

As- 7 819123 71 £H ' ^ ~~ ~ ES559 - V ^ ' 

where we omit couplings that are zero at this point, lies on a limit cycle. Among the zero 
couplings only the imaginary part of 2/5 and the real part of y% have nonzero beta functions 
and are thus generated on the limit cycle. Since not all imaginary parts of y\ 7 ... t s can be 
rotated away, the theory violates CP. For the entries of Q and P we find 

ttl 3 a/891 563478 -2 356 536^19 370, 

o {3) = — - 2 061664 + 143 986ci + 127268c 2 

H 3 763 549 370 814194 v 



- 735 868c 3 + 63 634c 4 - 735 868c 5 - 1 117968c 6 - 1 593 120c 7 

+ 654 696c 8 + 1 309 392c 9 + 1 726 320ci + 2 146 7B2c n - 25 316 928c 12 
+ 24 431 904ci 3 - 863 136c 14 + 4 779 648ci 5 + 106 491c 16 

- 212 982c 17 + 212 982c 18 + 106 491ci 9 - 212 982c 20 ), 



(2.5) 



and 
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(2.6) 



f3) 18\/297 187 826 - 785 512^19 370 , 

» = - 1881774685407097 (389632 + 43 °° C ' + ^ ~ W5U * C 

+ 25 360c 4 - 105 124c 5 - 94 632c 6 - 357 744c 7 + 93 528c 8 + 187 056c 9 

+ 276 648c 10 + 276 648c n - 3 616 704ci 2 + 3 490 272ci 3 - 155 844ci 4 

+ 862 992c 15 + 15 213c 16 - 30 426c 17 + 30 426ci 8 + 15 213ci 9 - 30 426c 20 ) - pf\ 

where the coefficients ci i ... j2 o are given by the contributions of the three-loop diagrams of 
Fig. [T] to the Yukawa beta function. None of the three- or four-loop contributions to (3 9 
appears in g^ 3 ) or pf\ For the other entries of P we find p^ 67s = 0, and that p 2 3 ] 4 are 
undetermined with p^ = — The condition for absence of anomalies of the dilatation 
current, TrP = 0, is thus p± +P3 = 0. We remark that and pf^ can be determined 
simply because the running couplings luiy 5 and Rey 8 run through zero at the point (2.4). 



4 Undctermined entries of P multiply operators that are conserved, i.e., they correspond to global sym- 
metries in the fermionic sector of the theory. 
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D 



(3) 
13 



D[ 4 (and its symmetric) 



D 



(3) 
15 




D (S) 




(and its symmetric) 




D\° g ' (and its symmetric) 




D 



(3) 
19 




L>2o (and its symmetric) 



Fig. 1: Three-loop diagrams that contribute to and p®. 



Note that both and pj[ receive contributions from exactly the same diagrams, 
although with different weights, and that twelve of these diagrams (D^-D^J , , and 
-Di3^) are exactly the diagrams that contributed to the frequency of the cycle of Ref. j3j. 
All diagrams involve at least two types of couplings, as expected from the "interference" 
arguments of Wallace and Zia 19 , as was also seen in our three-loop calculations in d — 4— e 
spacetime dimensions [3]. 

In dimensional regularization with d = 4 — e the three-loop diagrams of Fig. [l] have 
simple e-poles and thus they contribute to the Yukawa beta function. The residues of the 

(3) 

simple e-poles of D\_ 20 lead to the coefficients ci^.^o in 



(167T 2 ) a /3 a |y D Ci(2/ 6 y*y d 2/*y e ) y AaMe + • ■ • + c 20 g 2 [(y b y* c t* A t* A y d ) ij + {% o j}]X ab cd 



s 



as explained, for example, in 20 . We performed the three-loop computation with the 
method developed in Ref. 21 and the results of Ref. 22 . Since g^ 3 ^ and pf' are gauge- 



invariant, we can easily incorporate a quick check in our calculation by using the full gluon 
propagator, with the gauge parameter £. We fincQ 



C\ 


= 3, 


c 2 


= -1, 


c 3 


= 2, 


C 4 = 


5. 


c 5 = 


C7 


i 

2' 


c 8 


3 

2' 


c 9 


i 

2' 


ClO — 


5 

8 • 


Cll = 


Cl3 


1 

16' 


C14 


= 3, 


Cl5 


3 
8' 


Cl6 — 


-7 + 3£, 


Ci 7 = 



1 

2 ' 


c 6 = 


3 

2' 


5 




5 


8 ' 


Cl2 = 


32 


4£, 


CIS = 


-7-£, 



ci 9 = 19 + 5£, c 20 = 



Inserting these into the expressions (2.5) and (2.6) we obtain 



m 20 745^891 563 478 - 2 356 536 v / 19370 = __ fi 
y 99 040 772 916163 

and 

(3) (3) 

Pi = ~Ps ■ 

That g^ 3 -* ^ indicates that we have a limit cycle in the RG running of a four-dimensional 
unitary, renormalizable, well-defined gauge theory. This is the first example ever exhibited 
of such behavior. As expected, there is no ^-dependence in the final answer. As expected, 
the dilatation current is automatically non- anomalous. These are nontrivial checks and 
lend credibility to our calculation. We have found in the same theory a distinct second 
limit cycle, in another position in coupling space, with exactly the same properties as the 
one we presented above. 

We have verified that our results can be trusted in the e — y ^ limit. More specifically, 
the expansion parameters are bounded on the cycle: | A | /16vr 2 < 5%, \y\ 2 /16ir 2 < 1%, and 
g 2 /l6ir 2 = 0.46%. Hence, they remain perturbative along the whole cycle. 

The only unsatisfactory feature of our example is the fact that, as can be seen from 



Eqs. (2.4), the tree-level scalar potential is unbounded from below. Still the model can be 
studied in perturbation theory, since the vacuum state <ft — is perturbatively stable and 
its non-perturbative lifetime r is exponentially long, ln(r) ~ l/max(A Q ). This is similar in 
spirit to perturbative studies of renormalization for <f) 3 models in six dimensions. However, 
we expect that four-dimensional limit cycles with bounded scalar potential also exist. Our 
expectation is based on our results in d = 4 — e, where by progressing from the simplest 



5 The symmetry factors are included in the c's. Diagrams -Dg 3 ' n have symmetry factor s = h, diagrams 
13 have s = 4, and diagram has s = i. All other diagrams have s = 1. 
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examples, which displayed unbounded tree-level potentials, to more involved examples, we 
found limit cycles with bounded tree-level scalar potentials (lj|3). In any case, the behavior 
of the effective potential in any of these theories remains an open question. 



3. Conclusion 

The existence of limit cycles brings to light a new facet of unitary four-dimensional QFT. 
Many new questions arise: 

• What is the nature of RG flows away from limit cycles? Are there flows to or from 
fixed points from or to cycles or ergodic trajectories? 

• Are there limit cycles in supersymmetric theories? 

• Are there limit cycles in d = 3 and d > 4? Are there strongly-coupled limit cycles in 
d = 3 that correspond to the e — y 1 limit of the d = 4 — e perturbative models? 

• Are there limit cycles one can be establish in more indirect ways, i.e., without the 
need of three-loop computations? 

• Are there new possibilities for beyond the standard model physics associated with 
? 
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What is the holographic description of limit cycles? (This question has been consid- 
ered in Refs. |24|.) 



• Are there applications for condensed matter systems? 

Answers to these questions will allow a more complete understanding of QFT, and may 
lead to a new class of phenomena with unique characteristics. It should already be clear, 
though, that RG flows display behavior that is much richer than previously thought. 
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